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Abstract 

In the context of our group's recently developed emergent quantum mechanics, and, in particular, 
based on an assumed sub-quantum thermodynamics, Max Planck's assumption of energy quanti- 
zation is explained by means of a purely classical physics. Moreover, under the same premises, also 
the energy spectrum of the quantum mechanical harmonic oscillator is derived. 
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1. INTRODUCTION 

In references , the Schrodinger equation was derived in the context of modelling 
quantum systems via nonequilibrium thermodynamics, i.e., by the requirement that the 
dissipation function, or the time-averaged work over the system of interest, vanishes iden- 
tically. The "system of interest" is a "particle" in terms of a harmonic oscillator embedded 
in a thermal environment of non-zero average temperature (i.e., of the "vacuum"). In more 
recent papers, we have illustrated the "particle" more concretely by using the concept of 
a "bouncer" (or "walker", respectively) gleaned from beautiful experiments by Couder's 
group 3j-|6|. Thus one assumes that the thermal environment is oscillating itself, with the 
kinetic energy of these latter oscillations providing the energy necessary for the "particle" 
to maintain a constant energy, i.e., to remain in a nonequilibrium steady state. Regarding 
the respective ("zero point") oscillations of the vacuum, one simply assumes the particle 
oscillator embedded in an environment comprising a corresponding energy bath. 



2. A CLASSICAL OSCILLATOR DRIVEN BY ITS ENVIRONMENT'S ENERGY 
BATH: THE "BOUNCER" 

Let us start with the following Newtonian equation for a classical oscillator 

m'x = —iiiuJqX — 2'jmx + F cosut . (2-1) 

Eq. (12.1 p describes a forced oscillation of a mass m swinging around a center point along x(t) 
with amplitude A and damping factor, or friction, 7. If m could swing freely, its resonant 
angular frequency would be ujq. Due to the damping of the swinging particle there is a need 
for a locally independent (white noise) driving force Fq. 

We are only interested in the stationary solution of Eq. (12. ip . i.e., for t ^> 7 -1 , where 7 -1 
plays the role of a relaxation time, using the ansatz 

x{t) = Acos(ut + cp) . (2.2) 

After a short calculation we find for the phase shift between the forced oscillation and the 
forcing oscillation that 

tan<£ = ^ U , (2.3) 

OJq — u z 
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and for the amplitude of the forced oscillation 

A(u) = - Fo/m . (2.4) 

To analyse the energetic balance, we multiply Eq. (12.1 p with x and obtain 

mxx + rmo G xx = —2m^x 2 + F cos(ut)x , (2.5) 

and thus, 

— ( -mx 2 H — mu 2 x 2 J = ~2m^x 2 + F cos(wt)x = , (2.6) 

where the first term within the brackets is the kinetic energy and the second one is the 
potential energy of the oscillator system. Therefore, the time-derivative of the sum must be 
zero since the sum is the whole energy of the system and has hence to be of constant value. 

Due to the friction the oscillator looses energy to the energy bath, represented by ^rwyx 2 , 
whereas Fq cos(ojt)x represents the energy which has to be regained from the energy bath 
via the force F cos cut. The two terms on the right-hand side must also be zero and we can 
hence write down the net work-energy that is taken up from the bouncer in the form of heat 
during each period r as 

t+T t+T 
t t 

Inserting A from Eq. (12 .4p . we can see that dW/ du = at u = Uq, in which case the energy 
reaches its maximum. 

To derive the stationary frequency u, we use the right-hand side of Eq. ( 12. 6p together 
with Eq. (12.21) to first obtain 

2m r )x = —2m"/Au)sin(u)t + if) = F cos cut . (2-8) 

As all factors, except for the sinusoidal ones, are time independent, we have the necessary 
condition for the phase given by 

— sin(wt + ip) — cos cut =^ (f = — — + 2nir (2.9) 

for all n£Z. Substituting this into Eq. ( 12. 3h . we obtain 



tan + 2mr =±00 = ^ — -, (2.10) 

2 / Un—OJ 
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and therefore 

u = u . (2.11) 

Therefore, the system is stationary at the resonance frequency uq of the free undamped 
oscillator. With the notations 

r = — , r := A(u ) = , (2.12) 

ujq Z^muo 

we obtain 

Wbouncer = W houncer (u ) = "fmu$r 2 T = 27ijmu r 2 . (2.13) 

The Hamiltonian of the system is the term within the brackets of the left-hand side of 
Eq. ([23D, 

H = -mi 2 + -mojQX 2 = const. . (2-14) 

We introduce the angle 6{t) := u t and substitute Eq. ( 12. 2 p into Eq. ( I2.14p . thus yielding 
the two equations 

f - rO + u 2 r = , (2.15) 

and 

r9 + 2f9 = 0. (2.16) 

From Eq. ( I2.16p . an invariant quantity is obtained: it is the angular momentum, 

L(t) = mr 2 e\t) . (2.17) 

With 6{t) = uot, and thus 9 = cjq, the quantity of Eq. (12.17P becomes a time-invariant 
expression, which we denote as 

h := mr 2 u . (2.18) 

Thus, we rewrite our result (I2.13P as 

WWncer = . (2.19) 



3. BROWNIAN MOTION OF A PARTICLE: THE "WALKER" 

Now we concentrate on the motion of our "bouncer" within a more irregular environment. 
That is, we now assume that the "particle" is not only driven via harmonic oscillation of 
a wave-like environment, but that there is also a stochastic element in its movement, as, 
e.g., due to different fluctuating wave-like configurations in the environment. Therefore, 
our "particle's" motion will assume a Brownian-type character. The Brownian motion of 
a thus characterized particle (which we propose to call a "walker"), is then described by a 
Langevin's stochastic differential equation with velocity u = x and a force /(£), 

mil = —m(u + /(f) , (3.1) 

with a friction coefficient £■ The time-dependent force /(f) is stochastic, i.e., one has as 
usual for the time-averages 



/(f) = 0, f(t)f(t') = 4>(t-t'), (3.2) 

where 0(f) differs noticeably from zero only for t < Q^ 1 . The correlation time £ _1 denotes 
the time during which the fluctuations of the stochastic force remains correlated. 

The standard textbook solution for Eq. ( 13 . 1 j) in terms of the mean square displacements 
x 2 is given (e.g., in [7|) by the Ornstein-Uhlenbeck equation, 

£S=^(C|t|-l+e-<l«) . (3.3) 
Q z m 

Note that, on the one hand, for t <C and by expanding the exponential up to second 
order, Eq. (I3.3P provides that 

— t 2 = ^t 2 = u 2 t 2 . (3.4) 
mm 

On the other hand, for t ^> one obtains the familiar relation for Brownian motion, 



x 2 ~ 2Dt , (3.5) 

with the "diffusion constant" D given by 

D = ^. (3.6) 
Cm 
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To obtain a better understanding of Equations (13.51) and (13. 6p . we want to detail here 
how they come about. One usually introduces a coefficient A that measures the strength of 
the mean square deviation of the stochastic force, such that 

(j)(t) = \5{t) . (3.7) 

Since friction increases in proportion to the frequency of the stochastic collisions, there must 
exist a connection between A and (, respectively. One then solves the Langevin equation 
( 13.11) in order to find this connection. Solutions of this equation are well known from the 

Q 

Ornstein-Uhlenbeck theory of the Brownian motion [7J. 

Since the dependence of f(t) is known only statistically, one does not consider the average 

value of u(t), but instead that of its square, 

t t 

^) = e- 2 ^ [ dr [ dr'e^'V(T-T / )- + ^e- 2 ^ 

/ { m (3.8) 



2(m 2 v ' u 2(m 2 ' 

with uq being the initial value of the velocity. For t ^> the contribution of Uq becomes 
negligible, i.e., (~ l then plays the role of a relaxation time. We require that our particle 
attain thermal equilibrium (|l|, |2j) after long times so that due to the equipartition theorem 



the average value of the kinetic energy becomes, as usual, 

1 



-mu\t) = -kT . (3.9) 

Combining Eqs. (13 .8p and (13. 9ft . one obtains the Einstein relation 

A = 2(mkT . (3.10) 

Similarly, one obtains the mean square displacement of x(t) for t ^> £\ (Note that, 
even if we use the same character x as for the oscillating particle, the meaning is different: 
x(t) then signified a deterministic harmonic displacement of mass point m in the case of an 
oscillating particle ("bouncer"), whereas x(t) now means a stochastic random walk variable 
for the particle that carries out a Brownian motion ("walker").) Concerning the latter, one 
integrates twice to obtain the confirmation of our result (13.51) . i.e., 

t t 

*® = I dT I dT '2^ e ^ = ^ = 2Dt • (311) 

o o 
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with the diffusion constant turning out as identical to Eq. (j3.6p . i.e., 



A kT 

D = ^r^ = 7-- 3 - 12 



Now we remind ourselves that we have to do with a steady-state system. Just as with the 
friction ( there exists a flow of (kinetic) energy into the environment, there must therefore 
also exist a work-energy flow back into our system of interest. For its calculation, we multiply 
Eq. (13 . 1 p by u = x and obtain an energy-balance equation. It yields for the duration of time 
nr, with the natural number n > chosen so that nr is large enough to make all fluctuating 
contributions negligible, the net work-energy of the walker 

t+nr t+nT 

iy walker = J m(¥dt = m( J vHfj&t. (3.13) 
t t 

Inserting (13. 9p . we obtain 

kT 



Wwaikcr = nTm(u 2 (t) = riTm( — = nr(kT , (3.14) 

m 

where r = 2-k/uo equals the period of Eq. (12. 7p . which is chosen in order to make the result 
comparable with Eq. (I2.19p . The work-energy for the particle undergoing Brownian motion 
can thus be written as 

2tt 

^walker = U — (kT . (3.15) 

4. THE "WALKING BOUNCER": DERIVATION OF E = hu 

Let us summarize what we have achieved so far. We have for both systems, i.e., oscillator 
and particle in Brownian- type motion (or "bouncer" and "walker", respectively), obtained a 
net work-energy flow into each system, respectively, in order to compensate for the respective 
energy losses due to friction. There is a continuous flow from the bath to the oscillator, and 
vice versa. Moreover, and most importantly, during that flow, for long enough times nr, 
the friction of the bouncer can be assumed to be exactly identical with the friction of the 
walker. For this reason we directly compare the results of Eqs. (I2.19P and (I3.15p . 

^bouncer = ^walker , (4.1) 
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providing 



27T 

n2^h = n—(kT . (4.2) 

tu 



Now, one generally has that the total energy of a sinusoidal oscillator exactly equals 
twice its average kinetic energy. Moreover, despite having a nonequilibrium framework of 
our system, the fact that we deal with a steady state means that our oscillator is in local 
thermal equilibrium with its environment. As the average kinetic energy of the latter is 
always given, for each degree of freedom, by kT/2, one has for the corresponding total 
energy that .Etot = kT. Now, one can express that energy via Eq. (I4.2p in terms of the 
oscillator's frequency luq, and one obtains 

E tot = kT = ifwjQ . (4.3) 

7 

Assuming the same friction coefficient for both the bouncer and the walker, i.e., 7 = £, 
we obtain the energy balance between oscillator and its thermal environment as 

kT = hw , (4.4) 

with the total energy of our model for a quantum "particle", i.e., a driven steady-state 
oscillator system, being now derived as 

E tot = hw . (4.5) 

Moreover, if we compare Eq. (14. 4p with the Langevin equation (13. ip . we find the following 
confirmation of Eq. (I4.5p . First, we recall Boltzmann's relation between the heat applied 
io an oscillating system and a change in the action function 5S = 5 J E^ in dt, respectively, 
l|, l2j providing 

VQ = 2u V(5S) . (4.6) 
SS relates to the momentum fluctuation via 

V(6S) = 6p =: mu , (4.7) 

where, as usual, 

ft VP , s 

» = - 2 — ■ (4.8) 
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Thus one obtains with (14. 4p that the friction in Eq. (13. ip is given by 

h VP 

mCu = rri'-fu = —m / y — . (4.9) 

2m P 

Then, as P = P e~ s ®/ kT , with Eq. f !4.4j) for the thermal bath as above, and with Eq. (14. 6 p 
one obtains that 

C = 7 = 2u; , (4.10) 
and therefore the friction term becomes 

VP 

m(u = -hu)Q— = VQ . (4.11) 

Note that with Eqs. (14.41) and (I4.10p one obtains the expression for the diffusion constant 

„ kT h , . 

D = — = — } 4.12 
Qm 2m 

which is exactly the usual expression for D in the context of quantum mechanics. 

With Eq. (14. lip one can also introduce the recently proposed concept of an "entropic 
force" [8]. That is, with the total energy equaling a total work applied to the system, one 
can write (with S e denoting the entropy) 

2tt 



Ptot = 2P kin =: FAx = TAS C = J VQdr 

o 

AQ (circle) = 2 ( ^ 



, , , , ^0- (4.13) 

Eq. (I4.13P provides an "entropic" view of a harmonic oscillator in its thermal bath. First, 
the total energy of a simple harmonic oscillator is given as P to t — mr 2 bj\j2 =: ftwo/2. Now, 
the average kinetic energy of a harmonic oscillator is given by half of its total energy, i.e., by 



Pkin = w 2 Wg/4 = huo/4, which — because of the local equilibrium — is both the average 
kinetic energy of the bath and that of the "bouncer" particle. As the latter during one 
oscillation varies between and fkdo/2, one has the following entropic scenario. When it is 
minimal, the tendency towards maximal entropy will provide an entropic force equivalent to 
the absorption of the heat quantity AQ = huj^/A. Similarly, when it is maximal, the same 
tendency will now enforce that the heat AQ = tvjJo/A is given off again to the "thermostat" 
of the thermal bath. In sum, then, the total energy throughput E tot along a full circle will 



equal, according to Eq. (I4.13p . 2P kin (circle) = 2ftui /2 = hu . In other words, the formula 



E = ftojQ does not refer to a classical "object" oscillating with frequency Uq, but rather to a 
process of a "fleeting constancy" : due to entropic requirements, the energy exchange between 
bouncer and heat bath will constantly consist of absorbing and emitting heat quantities such 
that in sum the "total particle energy" emerges as Hcoq. 



5. ENERGY SPECTRUM OF THE HARMONIC OSCILLATOR FROM CLASSI- 
CAL PHYSICS 

A characteristic and natural feature of nonequilibrium steady-state systems is given by 
the requirement that the time integral of the so-called dissipation function Q t over full 
periods r vanish identically [1]. Assuming that our oscillator has a characteristic frequency 
Uo = 2tt/t, one defines the dissipation function w.r.t. the force in Eq. ( 12. ip over the integral 



Here, we assume a generalized driving force F to have a periodic component such that 
F(t) oc e luJ(,t . Then one generally has that 



o 



dF oc e iuJo{t+T) - e iWot , (5.2) 
and so the requirement Eq. (15. ip generally provides that 

uj dt = 2mr , for n = 1, 2, . . . . (5.3) 



n 

(Incidentally, this condition resolves the problem discussed by Wallstrom [9] about the 
single-valuedness of the quantum mechanical wave functions and eliminates possible contra- 
dictions arising from Nelson-type approaches to model quantum mechanics on a "particle 
centered" basis alone.) 

So, we are dealing with a situation where a "particle" oscillates with an angular frequency 
ido driven by the external force due to the surrounding ( "zero-point" ) fluctuation field, with 
a period r = 1 = For the type of oscillation we have assumed simple harmonic motion, 
or, equivalently [10], circular motion, and we generally have that the total ("zero-point") 



energy is 



1 2 2 ^0 /p. „x 

E = -mr ui Q = —— . (5.4) 
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Then, for slow, adiabatic changes during one period of oscillation, the action function over 
a cycle is an invariant, 

S = j> p ■ dr = ^- j> mu r ■ dr , (5.5) 

with r = r • f . This provides, in accordance with the corresponding standard relation for 
integrable conservative systems jl|, i.e., 

dS = , (5.6) 



that 

S = ^mr 2 u . (5.7) 

However, the external driving frequency and the particle's frequency uo, respectively, are 
not just in simple synchrony, since one has to take into account also the type of energy 
exchanges of the "particle" with its oscillating environment. Generally, there exists the 
possibility (within the same boundary condition, i.e., on the circle) of periods r n = ^ = 
with n — 1, 2, . . ., of adiabatical heat exchanges "disturbing" the simple particle oscillation 
as given by Eq. (15. 4ft . That is, while we have so far considered, via Eqs. (15. 5ft and (15. 6ft . a 
single, slow adiabatic change during an oscillation period, we now also admit the possibility 
of several (i.e., n) periodic heat exchanges during the same period, i.e., absorptions and 
emissions as in (I4.13p . The action integrals over full periods then more generally become 

T T T 

j> dS(r n ) := - J Sdt = J E tot dt = hju dt. (5.8) 



Thus, one can first recall the expressions (15 .4j) and (15. 7L respectively, to obtain for the case 
of "no additional periods" (n = 0) the basic "zero-point" scenario 

h — 1 

S = - , and E = U = - faj Q . (5.9) 

Secondly, however, using ( 15. 3p . one obtains for n = 1,2,... that 

dS(r n ) = Inixh = nh . (5.10) 
This provides a spectrum of n additional possible energy values, 

E(n) = nhwo , (5-11) 
11 



such that, together with Eq. (15. 9p . the total energy spectrum of the off-equilibrium steady- 
state harmonic oscillator becomes 

- ^ = E(n) + U = E(n) + E = ( n + ]- ] Huj , with n = 0, 1, 2, . . . . (5.12) 
at \ 2 J 

Note that to derive Eq. (I5.12p no Schrodinger or other quantum mechanical equation was 
used. Rather, it was sufficient to invoke Eq. (15.21) . without even specifying the exact expres- 
sion for F . 
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